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I. INTRODUCTION 



Preserving the coherence of quantum states and controlling their unitary evolution is one of the fundamental goals 
of Quantum Information Processing [0. When the system Hamiltonian is invariant under particle permutations, 
the exchange operator Eij interchanging particles i and j is a constant of the motion, and definite symmetry of a 
state will be conserved. Models of quantum computers based on identical bosons or fcrmions must of course respect 



2fl that active quantum error correcting codes 
1 fail for identical particles in the presence of 



this elementary requirement. It was pointed out in a recent paper | 
(QECCs) H designed to correct independent single-qubit errors, wi 
exchange errors. The reason is that exchange acts as a two-qubit error which has the same effect as a simultaneous 
bit flip on two different qubits. Of course, QECCs dealing explicitly with multiple-qubit errors are also available, 
so that exchange errors can readily be dealt with provided one accepts longer codewords than are needed to deal 
with single-qubit errors Q. For example, in Ref. a nine-qubit code is presented which can correct all single-qubit 
errors and all Pauli exchange errors. This is to be compared with the five-qubit "perfect" code which protects (only) 
against all single-qubit errors @. While the ninc-qubit code is longer than the "perfect" code, it is shorter than a 
code required to protect against all two-qubit errors. 

A different error model which has been considered by several authors is that in which qubits undergo collective, 
rather than independent errors [0-0] . The underlying physics of this model has a rich history: it dates back at 
least to Dicke's quantum optics work on superradiance of atoms coupled to a radiation field, where it arose in the 
consideration of systems confined to a region whose linear dimensions are small compared to the shortest wavelength of 
the field Jl(| . The model was later treated extensively by Agarwal in the context of spontaneous emission JllJ . It was 
only recently realized, however, that in the collective decoherence model there exist large decoherence-free subspaces 
(DFSs), which are "quiet" Hilbert subspaces in which no environmentally-induced errors occur at all [0J|]. Such 
subspaces offer a passive protection against decoherence. Collective decoherence is an assumption about the manner 
in which the environment couples to the system: instead of independent errors, as assumed in the active QECC 
approach, one assumes that errors are strongly correlated, in the sense that all qubits can be permuted without 
affecting the coupling between system and bath. This is clearly a very strong assumption, and it may not hold exactly 
in a realistic system-bath coupling scenario. To deal with this limitation, we have shown recently how DFSs can be 
stabilized in the presence of errors that perturb the exact permutation symmetry, by concatenating DFSs with QECCs 
|^| . Concatenation is a general technique that is useful for achieving fault tolerant quantum computation p^ , |i~3[ , and 
trades stability of quantum information for the price of longer codewords. It is our purpose here to analyze the effect 
of exchange errors on DFSs for collective decoherence. These errors are fundamentally different from those induced 
by the system-bath coupling, since they originate entirely from the internal system Hamiltonian. We will show that 
by use of the very same concatenation scheme as introduced in Ref. ^| (which was designed originally to deal with 
system-bath induced errors), a DFS can be stabilized in the presence of exchange errors as well. 

The structure of the paper is as follows. We begin by briefly recalling the origin of the exchange interaction in 
Sec. H and present some Hamiltonians modelling this interaction. We then present, in Sec. Ill, a short review of the 
Hamiltonian theory of DFSs. Next we discuss in Sec. |y| the simplest model, of constant exchange matrix elements, 
and show that DFSs are immune to exchange errors in this case. Our main result is then presented in Sec. when 
we analyze the effect of exchange errors in the case of arbitrary exchange matrix elements. We show that a DFS 
is invariant under such errors, and conclude that concatenation with a QECC can generally stabilize DFSs against 
exchange. 
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II. MODELLING EXCHANGE IN QUBIT ARRAYS 



The exchange interaction arises by virtue of permutation symmetry between identical particles, in addition to some 
interaction potential. Exchange is caused by the system Hamiltonian, and is unrelated to the coupling to an external 
environment. Exchange thus induces an extraneous unitary evolution on the system, but does not lead to decoherence. 
To model exchange it is sufficient to consider a Hamiltonian of the form 

1 K 

H CK = — Jij Eij , ( 1 ) 

where the sum is over all qubit pairs, Jjj are appropriate matrix elements, and 

Eij\ei, ...,€i, ...,e K ) = \ci, €j, £i, ex). (2) 

Eij thus written is a general exchange operator operating on qubits i and j of a -ftT-qubit state. 

Typical examples of Hamiltonian leading to exchange are |14f : (i) the Heisenberg interaction between spins 

■ffHcis = 2 X ' S 3 ( 3 ) 

where S; = (erf , erf, erf) is the Pauli matrix vector of spin i; (ii) the Coulomb interaction 

#Coui = g X) X) J i3 a L a ™' a ]*> a j° ( 4 ) 

i^i <r,<r' 

where a\ a (a,i a ) is the creation (annihilation) operator of an electron of spin a in Wannier orbital i. Jij is the exchange 
matrix element and is given for electrons by 

jC = _ g2 J dYdY , ™*{* ~ R>(r' ~ R^Vy - R>(r - Rj) ^ (5) 

where Ri is a lattice vector and w is a Wannier function . This is a rather generic form for the exchange matrix 
element; in other cases w would be replaced by the appropriate wave function and the Coulomb interaction e 2 /|r — r'| 
by the appropriate potential. The important point to notice is that the exchange integral depends on the overlap 
between the wave functions at locations i and j. Thus exchange effects generally decay rapidly as the distance |Rj — Rj | 
increases. An important simplification is possible when interactions beyond nearest neighbors can be neglected (i.e., 
Jij — if i and j are not nearest neighbors) in which case the approximation Jjj = J is often made. 

In the Coulomb case the interpretation of al cr ai er >aj a .,ajcr as an exchange operator is quite clear: spin a is destroyed 
at orbital j and is created at orbital i, while spin a' is destroyed at orbital i and is created at orbital j. The net effect 
is that spins a and a' are swapped between the electrons in orbitals i and j. In the Heisenberg case one can verify 
that the operator Si • also implements an exchange. Let / denote the identity operator, Xi the Pauli matrix af 
operating on qubit i, etc. A qubit state is written as usual as a superposition over a z eigenstates |0) and |1). Then, 
defining 

E^ = 1(1 + Si- Sj) = ^(I + Xi®Xj + Yi® Yj + Z t ® Zj) , (6) 
it is easily checked that Eq. (0) is satisfied [ ^5[ . 

III. REVIEW OF DECOHERENCE FREE SUBSPACES 

We briefly recall the Hamiltonian theory of DFSs J|,|l(| . Given is a system-bath interaction Hamiltonian 

H sb = J2f x ®B x , (7) 

A 

where Fx and Bx are, respectively, the system and bath operators. The decoherence free states are those, and only 
those states which are simultaneous degenerate eigenvectors of all system operators appearing in Hsb- 
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(8) 



The eigenvalues {c\} do not depend on The subspace spanned by these states is a DFS, meaning that under 
Hsb the evolution in this subspace is unitary, and there is no decoherence. This results in a passive protection 
against errors, to be contrasted with the active QECC approach. Of particular interest is the case where the {F\} 
are collective operators, such as the total spin operators 

K 

S a = ^af a = x,y,z. (9) 

i=l 

These operators satisfy su(2) commutation relations, just like the local af Pauli operators: 

[S a ,Sp] — 2ie a ji 1 S 1 . (10) 

This situation, referred to above as collective decoherence, arises when the bath couples in a permutation-invariant 
fashion to all qubits. In this paper we shall confine our attention to collective decoherence, and shall employ the 
term DFS exclusively in this context jL7). With a system-bath interaction of the form TJsb = ^2 a S a ® B a (as, 
e.g., in the Lamb-Dickc limit of the spin-boson model), a combinatorial calculation shows (see appendix) that the 

number of encoded qubits is log 2 K\j [{K/2 + 1)1(^/2)!] — 5° K — | log 2 K. The resulting decoherence free code 
thus asymptotically approaches unit efficiency (number of encoded qubits per physical qubits), and is therefore of 
significant interest. In the collective decoherence case, since the S a are the generators of the semisimple Lie algebra 
su(2), the DFS condition Eq. (||) is satisfied with c tt = |. This means that the decoherence free states {\j)} 
are su(2) singlets: they are states of zero total spin, and belong to the one-dimensional irreducible representation 
of su(2). For example, for K — 2 qubits undergoing collective decoherence, there is just one decoherence free state: 
(|01> - \1Q))/V2, i.e., the familiar singlet state of two spin 1/2 particles. For as few as K = 4 there are already two 
singlet states, spanning a full encoded decoherence free qubit f\. 



IV. DECOHERENCE FREE STATES AND EXCHANGE WITH CONSTANT MATRIX ELEMENTS 

A simple situation arises when we can assume that J^- = J/K for all i,j, i.e., without the restriction to nearest 
neighbor interactions. This long-range Ising model is thermodynamically equivalent to the mean-field theory of 
metallic ferromagnets, and there exist some examples of metals (e.g., HoRh4B4) that are well described by it ifLsf . 
At present the relevance of such materials to quantum computer architectures is not clear. We also stress that in the 
vast majority of physical examples exchange correlations decay exponentially fast with the distance between particles. 
The case of arbitrary exchange matrix elements is dealt with in the next section. We consider the long-range model 
here mainly for its simplicity and for the remarkable result that DFSs are completely immune to exchange errors in 
this case. 



We have for S = (S X ,S V ,S, 



S 2 = S ■ S = SKI + 2^Xi ® Xj + Y, L ® Yj + Z h (11) 



so that the exchange Hamiltonian can be rewritten as 

J K 

He * = 7^ zZ ( J + Xi X i + Yi ® y j + Zl ® Z ^ 

= ^[(K 2 -4K)I + S*]. (12) 

Whereas the DFS condition guarantees that no decoherence is caused by the coupling to the bath, uncontrolled 
unitary evolution due to the system Hamiltonian may still pose a significant problem. This is exactly the case in the 
presence of exchange errors, as described above. However, using Eq. ([ij]) and recalling that the DFS states have zero 
total spin, we see that in the collective decoherence case the DFS is in fact automatically protected against exchange 
errors: 



\rl>) = v\i/>), (13) 
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where \ip) is a DFS state and v = (J/K)(K 2 — 4fT)/8. Since the constant v does not depend on ip, this implies that 
under the unitary evolution generated by H ex , a DFS state accumulates an overall, global phase e wt . This phase 
is not measurable and does not affect the decoherence time. Thus in the Jy = J model a DFS does not undergo 
exchange errors, and the smallest DFS (K = 4 physical qubits) already suffices to encode a full logical qubit. 

V. DECOHERENCE FREE STATES AND ARBITRARY EXCHANGE MATRIX ELEMENTS 

We now analyze the effect of arbitrary exchange errors on DFS states for collective decoherence. We show that by 
concatenation with QECCs, DFSs can be stabilized against these errors. 

A. Decoherence Free Subspaces are Invariant Under Exchange 

The exchange operator commutes with the total spin operators. To see this, use the definitions of these operators 
in Eqs. (||) and (||), and let S l J = {^2k^i j a kj- Since they act on different qubits, S l J clearly commutes with £Iy. 
Now, using a a o^ = 8 a pl + ie Q /3 7 <7 7 : 

SaEy = [S a - (of + of)] By + (of + of) Ey 

\k=£i,3 / \ /3—x,y,z 

The last term in this expression vanishes since £ a/ 3 7 = —e ai p and we are summing over all /3, 7 values. Thus 

SaEii = SgEij + of + a? = E, tJ S a . (15) 

Now let be a decoherence free state (which it is for collective decoherence iff S a \ip) = ||). Since S a (Eij\ip)) = 
EijS a \ip) = 0, it follows that Eij\tp) also is decoherence free. We have thus proved: 

Theorem I. Let H be a decoherence free subspace against collective decoherence errors, and Eij an exchange 
operation on qubits i and j. Then EifH. = H. 

The significance of this result is that exchange errors act as errors on the encoded DFS qubits, i.e., they keep 
decoherence free states inside the DFS. The exact way in which these errors are manifested is a difficult problem. 
Exchange operations are transpositions in the language of the permutation group Sk , and are known to generate this 
group p9| . For a given number K of physical qubits the action of the exchange operators will realize a 2^-dimensional 
reducible representation of Sk ■ The DFS for collective decoherence on these K qubits is the set of one-dimensional 
irreducible subspaces in the irreducible representations (irreps) of Sk, which appear with multiplicity (k/2+i)\(k/2)\ 
(see appendix). For K = 4 the DFS is 2-dimensional (=multiplicity of the ID irreps), encoding one qubit. Therefore 
in this case exchange errors will act as the usual Pauli errors on a single (encoded-) qubit. Correction of exchange 
errors for K = 4 can then be done entirely within the DFS by using a quantum error correcting code for single-qubit 
errors. This observation naturally leads one to consider concatenating the DFS codewords with such a code, as done 
in the concatenated code of Ref. M. That paper showed that the concatenated DFS-QECC code can in fact deal with 
the more general case of both errors inside the DFS (as is our case here) , and errors that take states outside of the 
DFS. We investigate the correction of exchange errors in detail for the K — 4 case in the next subsection. For K > 4 
qubits, the dimension of the DFS is greater than 2 (e.g., for K = 6 it is 5), and the action of exchange errors will 
correspondingly be represented by higher dimensional irreps of Sk- To correct such unitary errors it will be necessary 
to resort to codes for "qufcits" (k > 2), such as stabilizer codes for higher-dimensional sytems ^0|, or polynomial 
codes |l3) . We defer the discussion of this case to a future publication |Q and focus here on the K = 4 case. 

B. Effect of Exchange Errors on the Four Qubit Decoherence Free Subspace 

Suppose that the qubits undergo collective decoherence in clusters of four identical particles, but different clusters 
are independent (as they might, e.g., in a polymer with an AAAABBBBAAAA... type of order). Each cluster 
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would then support a two-dimensional DFS, accommodating a single encoded DFS qubit. The K — 4 physical qubits 
DFS states can then be written as Jll| 

n,-wdi>. iD- g^iWHa , <i6) 

where 

\a) = |0110) + |1001), \b) = |1010) + |0101), |c> = |0011> + |1100>. (17) 

Note that the mutually orthogonal states \a), \b) and |c) are sums of complementary states. Moreover, the four 
qubits play a symmetrical role (i.e., and 1 appear equally in all four positions in both |0) and |1)). This dictates 
that exchange of qubits in symmetrical positions should have the same effect. In other words, we expect E12 to be 
indistinguishable from £34, and similarly for {£4.3, £24} and {£23, En} (although for a linear geometry most physical 
exchange mechanisms will yield | J23I > |«/i4|)- This expectation is born out; in the {|0), |1)} basis we find, using 
straightforward algebra, that the six exchange operators can be written as 



E12 = £ 



34 



-1 
1 



Evs = E 2A = R (tt/3) - -^-X + -Z 

AT _ 1 

E 14 = E 23 = R(-tt/3) = -^-X + -Z, (18) 

where R(0) = R(0)Z, and R(6) = f V Thus R(0) is a reflection about the x-axis followed by a 

counter-clockwise rotation in the x, y plane. In writing these expressions, the matrices operate on column vectors 
such that |6> = Q and |1) = (°), and X,Z are the encoded Pauli matrices, i.e., the Pauli matrices acting on the 
DFS states (and not on the physical qubits). Thus, exchange errors act as encoded Pauli errors on the DFS states. 

Using this observation, it is possible to protect DFS states against such errors by concatenation with a QECC 
designed to correct single qubit errors. The critical point is that this QECC will now correct single encoded qubit 
errors. This requires an additional encoding layer to be constructed. In particular, suppose we add such an encoding 
layer by using DFS qubits to build codewords of the five-qubit "perfect" QECC ||. These codewords have the form 
I £i) 1^2) 1^3) I £4) l^s}? where e = 0,1, and j in ij is now a cluster index. Since the five-qubit QECC can correct any 
single qubit error, in particular it can correct the specific errors of Eq. ( |l8| ) which the encoded DFS qubits would 
undergo under an exchange interaction on the physical qubits in a given cluster. However, the error detection and 
correction procedure must be carried out sufficiently fast so that exchange errors affecting multiple blocks at a time 
do not occur, or else concatenation with a code that can deal with t > 1 independent errors is needed. The typical 
time scale for exchange errors to occur is l/2\ Jij\, where Jij is the relevant exchange matrix element. 

This 20-qubit concatenated DFS-QECC code is precisely the one discussed in Ref. jjj, where it was shown that it 
offers protection against general collective decoherence symmetry breaking perturbations. Our present result shows 
that this concatenated code is stable against exchange errors as well. 

We note that it is certainly possible to find a shorter QECC than the five-qubit one to protect against the restricted 
set of errors in Eq. (|l^). However, such a code would not offer the full protection against general errors that is offered 
by concatenation with the perfect five-qubit code, and thus would not be as useful. 



VI. SUMMARY AND CONCLUSIONS 



To conclude, in this paper we considered the effect of unitary exchange errors between identical qubits on the 
protection of quantum information by decoherence free subspaces (DFSs) defined for a qubit array. We showed that 
in the important case of ideal collective decoherence (qubits are coupled symmetrically to the bath), for which a 
perfectly stable DFS is obtained, DFSs are additionally invariant to exchange errors. Thus such errors generate 
rotations inside the DFS, but do not take decoherence free states outside of the DFS. Consequently it is possible 
to use, without any modification, the concatenated DFS-QECC scheme of Ref. || in order to protect DFSs against 
exchange errors, while at the same time relaxing the constraint of ideal collective decoherence, and allowing for 
symmetry breaking perturbations. This is useful for quantum memory applications. Since exchange interactions 
preserve a DFS, an interesting further question is whether they can be used constructively in order to perform 
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controlled logic operations inside a DFS. We have found the answer to be positive, and that it is actually possible to 
perform universal computation in a fault tolerant manner inside a DFS for collective decoherence using only two-body 
exchange interactions p^] . 
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APPENDIX: DIMENSION OF DECOHERENCE FREE SUBSPACES FOR COLLECTIVE 

DECOHERENCE 



In view of the fact that the total spin operators 5 Q satisfy spin-1/2 commutation relations, it follows from the 
addition of angular momentum that the operators 5 2 and S z have simultaneous eigenstates given by 



S 2 \S,m) = S(S+l)\S,m), 



SzlS.m) = 



(Al) 



where m = ~S, -5 + 1, 5 and 5 = 0, 1, K/2 (for K even), 5 = 1/2, 3/2, K/2 (for K odd). The \S,m) states 
are known as Dicke states [^0|Jll| . The degeneracy of a state with given 5 is 



K\(2S+1) 



{K/2 + S+l)\(K/2- S)V 



(A2) 



which for 5 = 0, i.e., the singlet states, coincides with the dimension of the DFS for K qubits undergoing collective 
decoherence cited in the text. 

It is interesting to derive this formula from combinatorial arguments relating to the permutation group of K objects, 
which we will do for 5 = 0. The result follows straightforwardly from the Young diagram technique. As is well known 
(see, e.g., fy|), the singlet states of su{2) belong to the rectangular Young tableaux of K/2 columns and 2 rows. The 
multiplicity A of such states is the number of "standard tableaux" (tableaux containing an arrangement of numbers 
which increase from left to right in a row and from top to bottom in a column), which is also the dimension of the 
irreducible representation of the permutation group corresponding to the Young diagram riK/2,2 ( an empty tableau) 
of K/2 columns and 2 rows. This number is found using the "hook recipe" |fl9|| , where one writes the "hook length" 
gi (the sum of the number of positions to the right of box i, plus the number of positions below it, plus one) of each 
box i in the Young diagram: 



A"! 



1L=i 9i 



E.g., for ri c .2 the hook lengths are: 



and one finds, with c = K/2: 



c + 1 


c 


c-1 
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2 


c 


c-1 


c-2 




2 


1 



K\ 



2,2) 



(K/2 + l)\(K/2)V 



(A3) 



(A4) 



(A5) 



which is indeed the 5 = case of the general degeneracy formula, Eq. ( |A2| ) 
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